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", Inertial-range scaling behavior of high-order (up to order A'' = 51) structure functions of a passively 

■ advected vector field has been analyzed in the framework of the rapid-change model with strong 

' small-scale anisotropy with the aid of the renormalization group and the operator-product expansion. 

, It has been shown that in inertial range the leading terms of the structure functions are coordinate 

independent, but powerlike corrections appear with the same anomalous scaling exponents as for 
the passively advected scalar field. These exponents depend on anisotropy parameters in such a way 
, ^ that a specific hierarchy related to the degree of anisotropy is observed. Deviations from power-law 

. behavior like oscillations or logarithmic behavior in the corrections to structure functions have not 

' been found. 

' PACS numbers: PACS numbers: 47.10.+g, 47.27. -i, 05.10.Cc 
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I. INTRODUCTION 



Q 
U 

, Justification of the basic principles of the Kolmogorov-Obukhov (KO) phenomenological theory Q, |^ Q and the 
investigation of possible deviations from its conclusions within the framework of a microscopic model is one of the 
I ^ I main tasks in the theory of the fully developed turbulence and related models, e.g. stochastic magnetohydrodynamics 
(MHD). 

According to the KO theory, the following single-time structure functions in the inertial range (r^ r <C ?■/) 

> ; 

O. SN{r)^{[vr{K,t)-Vr{K\t)f), r=|x-x'| (1) 

'nT ■ 

, are independent of both the external (integral) scale r/ and internal (viscous) scale rd, the latter being tantamount to 
^ ' independence of viscosity. These requirements are the famous first and second hypotheses of Kolmogorov, respectively. 
, In Eq. ^ Vr denotes the component of the velocity field directed along the separation vector r = x — x'. Dimensional 
' arguments then determine the scale-invariant form of the structure functions ^ as 

.Si SN{r) ^ const X {erf , 

^ ' where e is the mean dissipation rate Q, H, Q . 

^ On the other hand, both theoretical and experimental results reveal some deviations from the KO theory 

viz. contradiction with the first Kolmogorov hypothesis. For the structure functions it means that they have to 
be modified in the following way 



SM{r) = {erf/'CN{r/rO, (2) 
where C,N{f/ri) are scaling functions with powerlike behavior in the asymptotic region r/r; ^ 1 



(Nir/ri) ~ const x j . (3) 

The singular dependence of the structure functions on r; in the limit r; ^ 00 together with nonlinearity of the 
exponents qn as functions of N is called "anomalous scaling". Theoretical explanation of such behavior is based on 
strongly developed fluctuations of the dissipation rate, i.e. intermittency Q,!^. 

A suitable and also powerful method of studying self-similar scaling behavior is that of the renormalization group 
(RG) 0, In ths theory of critical phenomena it successfully explains the origin of the critical scaling. This 

technique is also applicable to the theory of turbulence, see Refs. [H |M HI and references therein. 

The traditional approach to the description of fully developed turbulence is based on the stochastic Navier-Stokes 
equation j^. Within the RG method the second Kolmogorov hypothesis (independence of the viscous length r^) 
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was proved for structure functions for a variety of realistic random forces, (see Ref. [g)- This fact leads to the 
existence of the IR scaling (r :i> r^;) with definite "critical dimensions" 

^[Sn] = -J (4) 

with the customary convention A[r] = —1 for scaling dimensions. In the framework of the RG approach, critical 
dimensions like Q arise as coefficients in the differential RG equations. On the other hand, it is not possible to infer 
the form of the scaling functions CA'(^/n) in Eq. from the RG equations. 

The standard way to investigate the dependence of the scaling functions on the argument r/ri in the limit case 
r/ri is the utilization of the operator-product expansion (OPE), see Refs. 0, B 0- The OPE leads to the 
following representation of scaling functions Q 

where summation over all possible composite operators F (i.e. products of fields and their derivatives) is implied (for 
details, see Sec. IIV|I . Ap are their critical dimensions, and coefficients Cp are regular functions of r,~^. 

Contrary to the theory of critical phenomena, where there are physical reasons to believe that all relevant composite 
operators have positive critical dimensions {Ap > 0) for physical values of parameters [this is why the leading 
term in the expansion is given by the trivial operator F = 1 (Ai = 0) in that case], in the theory of fully 
developed turbulence based on the stochastic Navier-Stokes equation critical dimensions of many composite operators 
are definitely negative for physical values of parameters. Existence of these "dangerous operators" leads to singular 
behavior of structure functions in the limit r/r; ^ 8J. In the stochastic Navier-Stokes model dangerous operators 
enter into the OPE in the form of infinite families with the spectrum of critical dimensions unbounded from below, 
and a nontrivial problem of the summation of their contributions arises. This is an unsolved problem of the theory. 

In a situation where there are difficulties to study the anomalous scaling in the stochastic Navier-Stokes model 
(this applies to the stochastic MHD as well) it does not seem to be unreasonable to consider simpler models, which 
have features similar to real turbulent flow, and, on the other hand, are easier for investigation. An important role 
in this study was played by the model of passive advection of a scalar quantity (temperature or concentration of 
the tracer) by a uncorrelated- in-time Gaussian velocity field Models of passively advected vector fields are 
straightforward generalizations of the model of the passive advection of a scalar field. 

In the present work the spatial structure of correlations of fiuctuations of the magnetic (vector) field b in a given 
turbulent fluid in the framework of the kinematic MHD Kazantsev-Kraichnan model (KMHD) is studied. These 
fluctuations are generated stochastically by a Gaussian random emf and a white in time and anisotropic self-similar 
in space Gaussian drift. The main goal is the calculation of the anomalous exponents as functions of the anisotropy 
parameters of the drift. From the mathematical point of view the present model is found to be similar to the model 
of a passive scalar quantity advected by a Gaussian strongly anisotropic velocity field in that the fluctuation 
contributions to the critical dimensions Ap oi the OPE representation © of the structure functions coincide in both 
cases and thus the hierarchical dependence on the degree of anisotropy is also the same. Here, numerical calculation 
of the critical dimensions Ap in the one-loop approximation has been extended to dimensions related to structure 
functions of order A'^ = 51 to explore possible departures from powerlike asymptotic behavior. However, contrary 
to the scalar case, in the inertial range the leading terms of the structure functions of the magnetic field themselves 
are shown to be coordinate independent with powerlike corrections whose exponents are generated by the calculated 
critical dimensions. 



II. KINEMATIC MHD KAZANTSEV-KRAICHNAN MODEL 



Consider passive advection of a solenoidal magnetic field b = b(x, t) in the framework of the KMHD model described 
by the stochastic equation 

dth = i/oAb (v • V)b + (b • V)v + f , (6) 

where dt = d/dt, A = is the Laplace operator, i/q is the coefficient of the magnetic diffusivity, and v = v(x, i) 
is a random solenoidal (owing to the incompressibility) velocity field. Thus, both v and b are divergence-free vector 
fields: V-v==V-b = 0. A transverse Gaussian emf flux density f = f (x, t) with zero mean and the correlation 
function 



Dj^ = (/,(x,i)/,(x',t')) -5(t-t')C,,(r/L), r = x-x' 



(7) 
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is the source of the fluctuations of the magnetic field b. The parameter L represents an integral scale related to the 
stirring, and Cy is a function finite in the limit L ^ oo. In the present treatment its precise form is irrelevant, and 
with no loss of generality, we take Cij(O) = 1 in what follows. The random velocity field v obeys Gaussian statistics 
with zero mean and the correlation function 

DU.,t) . (.,(x,t).,(0,0)) d^^^^^^^^S0^,^ (8) 

where r; is another integral scale. In general, the scale r; may be different from the integral scale L, below we, 
however, take ri ~ L. Dq > is an amplitude factor related to the coupling constant go of the model by the relation 
Dq/vq = go = A'^, where A is the characteristic UV momentum scale, and < e < 2 is a free parameter. Its "physical" 
value e = 4/3 corresponds to the Kolmogorov scaling of the velocity correlation function in developed turbulence, d 
is the dimensionality of the coordinate space. In the isotropic case, the second-rank tensor Tij(k) in Eq. ||SJ| has the 
simple form of the ordinary transverse projector: Tij(k) = -Pij(k) = Sij — kikj/k"^ . 

In what follows we will be interested in the asymptotic behavior of the structure functions Sjq (r) within the inertial 
range [defined by the inequalities rd ^ r <^ ri, where Vd — is an internal (viscous) scale], which represent the 
equal-time correlations of the A^th power of the projection of the field b onto the direction along the separation vector 
of two different space points x and x' 



Dimensional analysis yields 



SM{r) = {[hrM-h.r{^,t)f), r=|x-x'|. (9) 



SN{r) = v~'''\''RN{r/rd, r/n), (10) 



where are functions of diniensionless parameters. When the random source field f and the velocity field v are 
uncorrelated, the odd functions S2n+i vanish, however. The standard perturbation expansion (series in go) is ill suited 
for calculation of structure functions 11U|) in the limit r/r^j — > oo and r/n — > 0, due to the singular behavior of the 
coefficients of the expansion. Therefore, to find the correct IR behavior it is necessary to sum the whole series. Such 
a summation can be carried out within the field-theoretic RG and OPE. A compact description of this procedure is 
presented in Refs. [l^ lisl Il6l| (see also Ref . [isj ) . Below we remind basic ideas and results referring to the isotropic 
case for simplicity of notation. 

The RG analysis can be divided into two main parts. First, the UV renormalization of structure functions ^ is 
carried out. As a consequence of this the asymptotic behavior of functions like © for r/vd 'S> 1 and arbitrary but 
fixed r/ri is given by IR stable fixed point (s) of the corresponding RG equations and for functions Q the following 
asymptotic form is obtained 

-N/2M 



where the scaling functions (,N{r/ri) remain unknown. In the standard language of the theory of critical phenomena 
(see for example Ref. the critical dimensions A[S'Ar] of the functions Sn are given by the relations A[S'Ar] = 

— where 7j!^ are "anomalous dimensions". The dimensions A[5jv] are calculated as series in e, and this is why 
the exponent e here plays the role analogous to the parameter 4 — d in the RG theory of critical phenomena. Another 
parallel is related to the parameter ri which is an analog of the correlation length Tc 0, 0] . 

Second, the small r/ri behavior of the functions ^N{r/ri) has to be estimated. This may be done using the OPE, 
which leads to the following asymptotic form in the limit r/ri — > 



where Cpir/ri) are coefficients regular in r/r;. The summation is implied over all possible renormalized scale-invariant 
composite operators F, and Ap are their critical dimensions. 

In the limit L/r — > 00 correlation function JTJ) of the random source field is uniform in space, which - as usual 
in stochastic models describing turbulence 0, IE ^1 0, E, E| - brings about composite operators with negative 
critical dimensions (dangerous composite operators) in the asymptotic analysis. Contributions of these dangerous 
operators to the OPE imply singular behavior of the scaling functions in the limit r/n — > 0. The leading term 
is given by the operator with the most negative critical dimension Ap. In our model the leading contributions to 
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even structure functions Sn are given by scalar operators = {bibi)^/'^ with their critical dimensions A^r, which 
eventually determine the asymptotic behavior of the structure functions S'^y of the form 



r, / ^ -N/2 N ( ^ 

SN[r) oci^o ' r — 



(11) 



where Ajv = —N + 7^^. Calculation shows that the anomalous dimensions 7]^ as well as the critical dimensions A^v 
in the model considered in the one-loop approximation are related to critical dimensions of composite operators of a 
simpler model of passively advected scalar field 0], viz. 7^ are given by 



7^ = ^ 



Ne 
2" 



N{N/2 - l)e 
{d + 2) 



(12) 



From relation Hll|l it follows that in the inertial range the structure functions are flat: S'Ar(r) ~ constl Taking 
into account the nonlinear dependence on TV of the anomalous dimensions H12|l we see that the deviation from the 
Kolmogorov scaling is extremely large in this case. It should be emphasized that this conclusion persists to all orders 
in perturbation theory and is not an artefact of the present one-loop calculation. Below it will be shown that these 
relations are stable against small-scale anisotropy. 

In the anisotropic case we will assume that the statistics of the velocity field is anisotropic at all scales and replace 
the ordinary transverse projection operator in Eq. JHJ) with the operator 



1 + aio- 



n • k 



Pij(k) + a2onsniPisi'k)Pji{'k) , 



(13) 



where Pij (k) is the usual transverse projection operator, the unit vector n determines the distinguished direction, 
and aio, a2o are parameters characterizing the anisotropy. The positive definiteness of the correlation function © 
imposes the following restrictions on their values: aio ,0120 > — 1- The operator (|13|l is a special case of the general 
transverse structure that possesses uniaxial anisotropy: 



T,,(k) = a(^)Py(k) + b{iP)nsniPMPji{k) , 



(14) 



where ■0 denotes the angle between the vectors n and k (n • k = fc cos ip) ■ Using Gegenbauer polynomials ^3 the 
scalar functions in representation (I14|l may be expressed in the form 



j=0 



aiP2iicosip) , b{il}) 



E 

1=0 



b.iP2i{cOSll)) . 



For the case of passively advected scalar, it was shown in Ref. that all main features of the general model with 
the anisotropy structure represented by Eq. are included in the simplified model with the special form of the 
transverse operator given by Eq. (|13|l . The same argument applies for the present case of passively advected vector 
field as well. 

The uniaxial anisotropy projector 1131) has already been widely used in analyzes of the anisotropically driven 
Navier-Stokes equation, MHD turbulence equations and passive advection equations However, these studies were 
limited to the investigation of the existence and stability of the fixed points with the subsequent calculation of the 
critical dimensions of the basic quantities leaving the calculation of the anomalous exponents in those models an open 
problem. 

The strong small-scale anisotropy (|14|l does not change the leading inertial-range term of the structure functions 
(|11|1 . but the anisotropy shows in the corrections to it. Indeed, combining the results of multiplicative renormalization 
and OPE in the manner sketched above for the isotropic case, we arrive at the conclusion that the inertial-range 
asymptotics of the structure functions of the passively advected vector field Sn is a constant independent of r with 
growing powerlike corrections effected by the small-scale anisotropy: 



SN{.r) 



N-Ne/2 / rd 



CN 



N / ^ ^ 

p=i ^ ' 



[«,p] 



A[M,p]+A[x,9] — Ajv 



where the critical dimensions Apf are given by Eq. H12|) at the order 0(e), whereas the critical dimensions A[jv,/ p] will 
be defined and calculated below [Eq. (|40|l ] . This is similar to the asymptotic behavior of the structure functions of 
the passive scalar advected by a compressible vector field jlSj] . 
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III. FIELD-THEORETIC FORMULATION, RENORMALIZATION, AND RG ANALYSIS 

The stochastic problem ©-(IHll is equivalent to the field-theoretic model of the set of the three fields $ = {b', b, v} 
with the action functional 

= i h'DbW + h'[-dt - (v V) + iyaA]h + b'(b ■ V)v - ^vl^-^v, (15) 

where b' is an auxiliary field (all required integrations over space-time coordinates and summations over the vector 
indices are implied). The first five terms in Eq. (I15|l represent the De Dominicis-Janssen action corresponding to the 
stochastic problem at fixed v (see, e.g., Refs. 22]), whereas the last term represents the Gaussian averaging over v. 
Z?fc and Dy are the correlation functions Q and ||SJ), respectively. 

In this field-theoretic language, the structure functions © are defined as 

SNir) = J P$[6,(x, t) - 6,(x', i)]^e^(*) (16) 

with the action S'($) defined above. 

Action IjlSI) is given in a form convenient for application of the quantum-field perturbation analysis with the standard 
Feynman-diagram technique. The quadratic part of the action determines the matrix of bare propagators. For the 
fields b' and b the propagators in the wave-vector-frequency representation are 



Cy(k) 



where Cij(k) is the Fourier transform of the function Cij{r/L) from Eq. 0. The bare propagator of the velocity 
field (vv}o = (vv) is defined by Eq. © with the transverse projector given by Eq. (|13|l . The interaction in the 
model is given by the nonlinear terms — ^((v • V)bi + 6J;(b • V)ui = b[VijiVjbi with the vertex factor which in the 
wave-number-frequency representation has the following form 

Viji = i{6ijki - Siikj). 

With the use of the standard power counting ^ (see also Ref. jTsj for peculiarities of rapid-change passive advection 
models) correlation functions with superficial UV divergences may be identified. These are correlation functions 
containing frequency- wave- vector integrals divergent in the limit e ^ with divergences brought about by integration 
over large wave numbers and correspondingly having non-negative wave-number dimension. In the present model 
superficial divergences exist only in the one-particle-irreducible (IPI) Green function Fb'b- In the isotropic case this 
Green function gives rise only to the renormalization of the term i^ob'Ab of action (|15(l and the corresponding UV 
divergences may be fully absorbed in the proper redefinition of the existing parameters go, so that all correlation 
functions calculated in terms of the renormalized parameters g and v are UV finite. 

When anisotropy is introduced, however, the situation becomes more complicated, because the IPI Green function 
Fb'b produces divergences corresponding to the structure b'(nV)^b in the action of the model [due to peculiarities of 
the rapid-change models the term (b'n)A(bn) possible on dimensional and symmetry grounds does not appear]. 
The term b'(nV)^b is not present in the original unrenormalized action H15|l . but has to be added to the renormalized 
action, therefore the model is not multiplicatively renormalizable. In such a case it is customary to extend the original 
action (|15|l by including all terms needed for the renormalization of the correlation functions and thus adding new 
parameters. As a result the extended model is described by a new action of the form: 

5($) = h,'Di,h' + h'[-dt - (v • V) + i^oA + Xo'^o(n • V)2]b + b'(b • V)v - ^vi^-^v, (18) 

where a new unrenormalized parameter xo has been introduced. 

Of course, the bare propagators H17(l of the isotropic model are modified and for the extended action (|18|l assume 
the form 

(6.5;.)o - (6;5,)S = . ^ , 1.^2 ' (19) 

■' ■' -iti; + j/ofc-^ + Xot'o(n • k)^ 



(^^^^■)o = I ■ ^ ( 1.^212 ' (20) 

] - z^ofc^ + Xot'o(n • k)^]^ 




FIG. 1: The (exact) graphical expression for the self-energy operator 'Etn, of the response function of the passive vector field. 
The plain line denotes the bare propagator H2()|l . and the line with slash (denoting the end corresponding to the arguments of 
the field b') corresponds to the bare propagator 11911 . 

After this modification all terms needed to remove the divergences are present in action (|18ll . therefore the model 
becomes multiplicatively renormalizable allowing for the standard RG analysis. The corresponding renormalized 
action may be written down immediately: 

5i?($) = ib'Ab' + h'[-dt - (v • V) + z/ZiA + xz/Z2(n • V)^]h + b'(b • V)v - ^vi^-^v . (21) 

Here, Zi and Z2 are the renormalization constants in which the UV divergent parts of the IPI response function Fb'b 
are absorbed. The renormalized action ()21|l leads to the multiplicative renormalization of the parameters voydo and 
Xo- 

1^0 = 1^2^, go = giJ.''Zg, xo = X^x' 

where v, g, and x a-re renormalized counterparts of the bare parameters, and /i is a scale-setting parameter with the 
same canonical dimension as the wave number. The anisotropy parameters aoi and aa2 are not renormalized, therefore 
their renormalized counterparts ai and a2 may be put equal to the unrenormalized parameters. In what follows, 
we will work in the minimal subtraction (MS) scheme, in which - in the one-loop aproximation - renormalization 
constants have the form 1 -I- A/e, where the amplitude A is a function of g, x, ai, 012, and d, but independent of e. 

Identification of the unrenormalized action ()18|l with the renormalized one H21|l leads to the following relations 
between the renormalization constants: 

Zl = Zu, Z2 = Zy^Zy^ Zg = Z^j^ . 

It has to be mentioned that the rapid-change models like H18|l have the nice feature that in all multiloop diagrams 
of the self-energy operator Sb'b closed circuits of the retarded bare propagators (bb')o are produced, because the 
propagator (vv)o is proportional to the 5 function in time. As a result, the one-loop self energy operator !]{,/{, with 
the graphical notation of Fig. ^is exact. 
The divergent part of the graph in Fig. ^ is 

Sb'bb) = -^^g^^{[(d-l)(d + 2)+ai(d+l) + a2]p2_(2ai-(rf2_2)a2)(n-p)2} , (22) 

where Cd = Sdli^ixf and Sd = 27r''/7r(d/2) is the area of the d-dimensional sphere of unit radius. The expression 
(|22|1 leads to a straightforward determination of the renormalization constants Z\ and Z^'- 

Z. = 1 - .^li'l-LV - IX"* + 2) + oi(<i + 1) + "2I ■ 

With the proper choice of the renormalization constants this renormalization procedure gives rise to UV-finite cor- 
relation functions with separate space and time arguments. Independence of the original unrenormalized model of 
the scale-setting parameter \x of the renormalized model yields the RG differential equations for the renormalized 
correlation functions of the fields, e.g. 

(P^ + ligdg -f I5^d^ - 7.^?.)(b(x, i)b(x', t))H - (24) 

with the definition T)^ = xdx for any variable x, and with the following definition of the RG functions (the j3 functions 
and the anomalous dimensions 7) 

7i = In Zi 

for any renormalization constant Zi , and 

I3g EE V^g = g{-e + 71) , (3-^ = V^x = x(7i - I2) , 
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where V = ^d^ denotes derivative with fixed bare parameters of the extended action (|18() . 
For the anomalous dimensions 71 and 72 we obtain h'om Eq. H23|l 

71 = 5jr^ P - !)((! + 2) + a,(<i + 1) + aj , (25) 

It should be emphasized that both the renormalization constants \2'6\ and the corresponding anomalous dimensions 
(|25|l and H26I) in the present model are exact, i.e., they have no corrections of order or higher. 
The fixed points X*) of RG equations are defined by the system of two equations 

/3g(5*,X*) = 0, /3x(5*,X*) = 0- 
The IR stability of a fixed point is determined by the condition that real parts of the eigenvalues of the matrix 



g=g* 
x=x* 



are positive. Calculation shows that the RG equations have only one non-trivial IR stable fixed point defined by 
expressions 

^* Cd[id-l){d + 2)+ai{d+l) + a2]' ^ ' 

-2a^ + (d^ - 2)a^ 

^* (d-l)(d + 2) + ai(d+l) + a2' ^ ' 

Both eigenvalues of the stability matrix lo are equal to e at this fixed point, therefore, the IR fixed point (|27|) . (|28|l is 
stable for e > and all values of the anisotropy parameters ai and a2- 

Rather unexpectedly, the P functions and, consequently, the fixed points of the present model of passively advected 
vector field are exactly the same as in the model of passively advected scalar field '13'|. In Sec. II VI it will be shown that 
this similarity is extended to the anomalous scaling dimensions of the composite operators in the OPE representation 
of the structure functions as well. 

The fixed point H27|l . H28(l governs the behavior of solutions of Eqs. H24|l and the like, and at large scales far from 
viscous length r ^ at any fixed ratio r/ri yields the scaling form 

(b(x, t)b(x - r, t)) = D^^^-'^Tln) , (29) 

for the unrenormalized correlation function (we remind that due to the absence of field renormalization renormalized 
and unrenormalized correlation functions are equal but expressed in terms of different variables). It should be noted, 
however, that the scaling function £,2{'''/ri) in Eq. is not determined by the RG Eqs. (|^ . 

This, however, is not enough to find the asymptotic scaling behavior of the structure functions because they are 
linear combinations of correlation functions with coinciding arguments. The latter contain UV divergences additional 
to those included in the renormalization constants H23|l . These additional divergences due to composite operators 
(products of fields and their derivatives with coinciding space and time arguments) can be dealt with in a manner 
similar to that applied to the divergences in the usual correlation functions . 



IV. RENORMALIZATION AND CRITICAL DIMENSIONS OF COMPOSITE OPERATORS 



A composite operator is any product of fields and their derivatives at a single space-time point x = (x, t), e.g. 
[b^(a;)]^ and [dihj{x)dibj{x)] . Structure functions Hlt)|) contain products of composite operators at two separate 
space points. These composite operators are integer powers of the field b^. Usually, the aim of renormalization of 
composite operators is to make UV finite all IPI correlation functions with the insertion of a composite operator F, 
i.e. quantities of the form 

(-F(/)ii(xi,ti) • • • (j)i^{Xm, tm))lPl ■ 
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Correlation functions with such insertions contain additional UV divergences, which also may be removed by a suitable 
renormalization procedure 0, ^] . Composite operators mix in renormalization, therefore an UV finite renormalized 
operator has the form = F + AF, where the counterterms ^-re a linear combination of the operator 

F itself and, in general, other unrcnormalized operators required to make all correlation functions - generated by 
the renormalized action - with the insertion of F^ UV finite. In this case homogeneous RG equations of the form 
()24|l may be obtained for certain linear combinations of renormalized correlation functions with composite-operator 
insertions. Such linear combinations (basis operators see below) exhibit IR scaling with definite critical dimensions 
AjF, whereas an arbitrary renormalized composite operator may be expressed as a linear combination of these basis 
operators. 

The general procedure is the following |EIE0)B0i0|: If {Pa} is a closed set of composite operators (i.e., they 
are mixed only with each other in renormalization), then the sets of renormalized and unrcnormalized operators are 
related through the matrix transformation: 

Fa = ^ ZapFp , (30) 

P 

7F = Zp^Vf^Zp , (31) 

where Zp = {Zap} is the renormalization matrix and 7f = {lap} is the corresponding matrix of anomalous dimensions 
for this set of operators. The renormalized composite operators obey the RG differential equations 

P 

which give rise to the matrix of critical dimensions Ap = {Aap} of the form 

Ap = d%- Atd^p +-i*p, At = -2 + e , (32) 

where dp and d'^, respectively, are the diagonal matrices of canonical wave- number and frequency dimensions of the 
operators (where the diagonal elements are sums of the corresponding dimensions of the operators included in the 
composite operator) and "^p is the matrix of anomalous dimensions at the fixed point H27|l . (jJHJ- 

Critical dimensions of the set F = {F^} are given by the eigenvalues of the matrix Ap. The basis operators pos- 
sessing definite critical dimensions are related to the renormalized composite operators by the matrix transformation 



P 

where the matrix Up = {Uap} is such that the transformed matrix of critical dimensions Ajr — UpApUp^ is diagonal. 

The structure functions contain, however, quantities which correspond to insertions of two composite operators. 
Therefore, it would seem that we would have to consider renormalization of products of two composite operators as 
well, the aim being then to render UV finite all IPI correlation functions with two insertions of composite operators. 
Superficially divergent correlation functions with operator insertions are identified by power counting similar to that of 
the basic renormalization. In the present model such a power counting shows that insertion of products of composite 
operators of the structure b'"(x, t)b"(x', t) does not bring about any new superficial divergences and it is thus sufficient 
to renormalize the composite operators themselves only in order to make the structure functions UV finite. Therefore, 
from the RG analysis of composite operators it follows - by virtue of relations (|30|l and (|33|l - that the structure 
function Sn may be expressed as a functional average of a quadratic form of basis operators: 

SN{r) = g (^x + ^r, < j (^x - ^r, ?j ^ (34) 

with coefficients Bap independent of spatial coordinates. Each term in expression (|34|l obeys the following asymptotic 
form in the limit rd 'ti r, r < ri 

(x+ ir,t) (x - lr,t))^ ^ Dt^^'K--^--r^rf^^h^, (^) (35) 

with the scaling functions still to be determined. 
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The physically interesting range of scales, however, is the inertial range, specified by the inequalities <C r <C r/. 
The limit r <C may be explored with the use of the OPE j5i,||^ as was already discussed in Sec. 13 The basic statement 
of the OPE theory is that the equal-time product of two renormalized composite operators can be represented in the 
form 



fS ( X + t\ (x - ir, < j = ^ C^p,{v)F^(^, t), (36) 
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where the functions Cafij are the Wilson coefficients regular in 1/ri, and are renormalized local composite operators 
which appear in the formal Taylor expansion with respect to r together with all operators that mix with them in 
renormalization. If these operators have additional vector indices, they are contracted with the corresponding indices 
of the coefficients Cap^- 

Without loss of generality we may take the expansion on the right-hand side of Eq. H36|l in terms of the basis 
operators with definite critical dimensions Ajf- The renormalized correlation function {F^'Fp')ij is obtained by 
averaging Eq. I|36|) with the weight generated by the renormalized action, the quantities appear now only on 

the right-hand side. Their asymptotic behavior for r/ri — > is found from the corresponding RG equations and is 
of the form (JF) oc ^ ■ Comparison of the expression for a given function {F^Ff)^ in terms of the IR scaling 
representation of correlation functions of the basis operators (I35|l on one hand and the OPE representation brought 
about by relation (|36|l on the other in the limit r; ^ oo allows to find the asymptotic form of the scaling functions 
Sq,^ (r/rj) in relation H35|l . 

The composite operators appearing in the expression for the structure function Stq are products of integer powers 
of the field hr of the form 6^~'"(x, t)6™(x', t). Thus, at the leading order in r their OPE contains operators of the 
closed set generated by the operator (x, t) . Power counting and analysis of the structure of graphs shows that 
this set of composite operators contains only operators consisting of exactly N components of the vector field b (i.e. 
no derivatives of the field components appear). Extracting the common scaling factor prescribed by the canonical 
dimensions l|32|) of these operators the basis-operator decomposition of any term of the structure function Stq may be 
written as 

V^"" (x+ir,t) C (x-ir,t)) X: A^,,{rln) (^)"""' (^)^' (37) 

where the coefficients Aaji-^{rlri) are regular in {r/riY. 

The decomposition H37|l reveals the inertial-range scaling form of the structure functions. The leading singular 
contribution in the limit r; — > oo, rd ^ is given by the basis operator T-y with the minimal critical dimension 
and operators J^a and with the minimal sum of anomalous dimensions + 7^- As a result, Sn have singular 
power-like behavior as r/ri — > 0: 

-In / ^ \ 



with the most negative exponent Aat in the basis set generated by the composite operator (x, t) and the most 
negative sum of anomalous dimensions 7j*jy^ -I- 7j*^ subject to the conditions M + K = N , p < M and q < K. 

In our case, the leading contribution to the sum H37|l from the OPE (|36|l will be given by the tensor composite 
operators constructed solely of the fields b without derivatives: bi^...hi {hihi)K It is useful to deal with the scalar 
operators obtained by contracting the tensor with the appropriate number of the anisotropy vectors n: 

F[iV,p](x,i) = [n.b(x,i)f [6.(x,t)6,(x,i)]' (38) 

with N = 2l+p. Power counting and analysis of graphs show that composite operators (|38|l for given N can be mixed 
only with each other in renormalization 

N 
p'=0 

therefore, the corresponding renormalization matrix Z^^ p][N' p'] is in fact block-diagonal, i.e., -^[at j,][Ar' p'j = for 
N'^N. 

A detailed account of practical calculation of the matrix of the renormalization constants Zi^N.p][N.p'] (which may 
be readily extended to investigation of all related problems) has been given in Ref. (13i] for the advection of a passive 
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scalar, therefore we will not describe all details of the determination of renormalization constants in the present vector 
model, rather we will discuss its specific features. 

It turned out that not only the (i functions in the vector and scalar models coincide, but the one-loop renormalization 
matrices as well. This nontrivial fact stems from the similarities of the mathematical structure of both models. In the 
model of scalar advection jl^i] the composite operators di^6...di 6{di9di9y constructed solely of the scalar gradients of 
the scalar admixture 9 are needed for calculation of the asymptotic behavior of the structure functions, whereas in our 
vector case the main contribution is given by composite operators constructed solely of the fields b without derivatives. 
As direct inspection of the relevant diagrams shows, the tensor structures arising upon functional averaging in both 
cases are in fact identical, which yields the same renormalization matrix Z[jv,p] [Ar,p'] in both models. Thus, it is not 
necessary to carry out complete calculations here. 

However, in Ref. 13] in the expressions for the general elements of the renormalization matrix of the composite 
operators there are misprints (for instance, in the definition of the quantity Qi in Eq. (76) of Ref. 13] Hi^ — Hq 
should be replaced by — Hq), although the numerical investigation of the critical dimensions is correct. Therefore, 
we present here the full formulae, in a slightly different form, however. 

The only nonzero elements of the matrix Z^N^p^^^y^ are 

__gCd_i _ 1 , gCd 1 ^ 

^[N,p][N,p-2] - ^2 _ 1 4^^1 ' Z[N,p][N,p] - 1 + ^2 _l 4^y2 , 

_ gC'd 1 gCd 1 

ZlN,p][N,p+2] - ^2 _ 1 4^^^ , ^[JV,p][JV,p+4] - ^2 _ 1 4^^" ' 

with the coefficients Qi defined as follows 

Qi = p{p -l)id+ l)(iJo(l + "2) + (Hi - H2){1 - 2ai + Saj) + H^^ai - as)) 
Q2 = Ho{{d - 2){N ~ pf - (1 + a2){d +l){p- l)p 

+ {N - p){i + d^ + 2d{p -l) + 2p + a2{d+ 1)(1 + 2p))) 
+ H2{{5 + 3a2 + ai{d - 2) - d){N - pf 

+ {d+l){\-ai+d + a2{d + 2)){p-l)p-{N -p){9-2d + d^ + 8p{d + 1) 

- ai(3 + + 2d{p - 1) + 2p) + 2a2(4 + d + hp{d + 1)))) 
+ iJ4(-(3 + 6a2 + ai{d - 5))(iV - pf + {d+ \){ai{d + 1) - a2{2d + 1) - d){p ~ \)p 

+ (A^ - p)(6(l + p{d + 1)) + 02(13 + d + 14p(d + 1)) 

-ai(9 + d2+8p + 2d(4p-l)))) 
+ H^{oi2 + ai)(3(iV - pf + d{d + 1)(?J - \)p - 6(iV - p){\^p{d + 1))) , 
Q3 = iJo(n - p)(-9 ~d^ + -Ip- d{-2 + N + p)+a2{d+l){n- 3(1 + p))) 
+ H2{N -p){d^ + d^{-^2 + N +p) + 6(6 - 3A^ + Ap) + (i(9 ~bN + lip) 

+ ^2(27 - 13A^ + 2lp + d^{l + 2p) + d(16 -IN + lip)) 

+ ai(9 + d2 -5A^ + 7p + rf(-2 + iV+j5))) 
+ Hi{N - p){-2{2 + d)(6 - 3A + (4 + d)p) 

+ ai{d^ +d^{-2 + N +p)+ 6(6 - 3A + Ap) + rf(9 - 5iV + IZp)) 

+ ^2(48 - 247V + Mp + d^(l + Ap) + d(25 - 12A + 26p))) 
+ He2{N - p){a2 - ai){2 + d){6-3N+{A + d)p) , 
Qi = Ho{3 + a2{d+ 1))(2 -N + p){N-p) 

+ H2{'iai - 6(2 + d) - 02(9 + 7d + d^)){2 - N + p){N ~ p) 
+ Hi(2 + d){-6ai + (1 + 2a2)(4 + d)){2 - N + p){N - p) 
+ He{ai - a2){8 + 6d + d^){2 -N+p){N-p), 
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where Hi are the functions 



Ho=2Fi(l,l/2,d/2,-x), 

i/2 =2i^i(l,3/2,d/2 + l,-x)/d, 

Ha - 2^^i(l, 5/2, d/2 + 2, -x):jJ-, 




Here, 2F1 is the Gauss hypergeometric function. The anomalous dimensions 1[n .p][N .p'] sre 

_ gCd ^ _ gCd 




Q 



'2 , 



'4, 



(39) 



and the matrix of critical dimensions H32|l is 



^[N,p][N,p'] — 




'^pp' + 7[Ar,p] [N,p'] J 



(40) 



where the asterisk stands for the value at the fixed point (|27|l . H28(l . This represents the critical dimensions of the 
composite operators (I38II at the first order in e. It should to be stressed that in contrast to the value of the fixed point 
(|27|l . (|28|l . which has no higher order corrections, the expressions for anomalous dimensions H39() have nonvanishing 
corrections of order and higher. 

The critical dimensions are given by the eigenvalues of the matrix (|40|l . As was already discussed in Ref. in 
the limiting isotropic case (ai = a2 = 0) this matrix becomes triangular, i.e., the eigenvalues are simply equal to the 
diagonal elements A[A^,p] = A[7v,p][Af,p]- 

Since our result for the anomalous dimensions is the same as in Ref. jlj] for the admixture of a passive scalar, all 
conclusions about the hierarchical behavior of the critical dimensions of the composite operators are also valid in the 
analysis of the present model. Nevertheless, the inertial-range asymptotic behavior of the structure functions in these 
two problems is completely different, because, first, in the scalar problem single-point products of the scalar are not 
renormalized, while in the vector problem they are, and, second, the leading contribution to the OPE is given by the 
products of derivatives of the scalar, whereas in the vector problem products of the field components themselves yield 
the leading contribution. 

In Ref. 13] the behavior of the critical dimensions A[Af, p] for N — 2,3,4,5, and 6 was numerically studied. The 
main conclusion is that the dimensions Ajv remain negative in anisotropic case and decrease monotonically as N 
increases for odd and even values of N separately. 

In present paper we concentrate our attention on the investigation of the composite operators (|38|l for relatively 
large values of N, namely we will analyze cases with N = 10, 11,20,21,30,31,40,41,50, and 51. Our aim has been 
to find out whether hierarchies which hold for small values of TV remain valid for significantly larger values of N, and 
the answer turned out to be in the affirmative. 

In Ref. ^3 several hypothetically possible structures of the matrix of critical dimensions (|40|l were discussed. In 
particular, the possibility that the matrix (|40|l for some ai and 02 would have a pair of complex conjugate eigenvalues 
A = Re A ± ilmA cannot be excluded a priori. In this case, the small-scale behavior of the scaling functions would 
have oscillating terms of the form 



with some real constants Ci , C2 . 

Another, in general, conceivable structure of the matrix H40|l is related to the situation when it cannot be diagonalized 
but only reduced to the Jordan normal form. In this case, the corresponding contribution to the scaling function 
would involve a logarithmic correction to the power-like behavior, viz. 




Rc A 



{Ci cos[(ImA)r/n] + C2 sin [(Im A)r/n]} , 




where A is the eigenvalue related to the Jordan cell. 
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FIG. 2: Behavior of the critical dimension A[10,p]/e for space dimension d = 3 and for representative values of p as functions 
of anisotropy parameters qi and a2- 




FIG. 3: Behavior of the critical dimension A[ll,p]/e for space dimension d = 3 and for representative values of p as functions 
of anisotropy parameters qi and 02- 



In Figs. I2I7I behavior of the eigenvalues of the matrix of critical dimensions A[Af,p] for relatively large values of the 
N are shown. It can be seen that only real eigenvalues exist in all cases, and also their hierarchical behavior discussed 
in Ref. 'l3'| is conserved. At first sight the curves for p = and p = 2 in the even case and the curves for p = 1 and 
p = 3 in the odd case in Figs. I4I7I appear to be crossing at the point ai = Q!2 = but in fact the curves are only 
visually running very near together at that point which is a mathematical consequence of the formulas for critical 
dimensions in the infinitesimal limit ai — > and a2 — > 0. 
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FIG. 4: Behavior of the critical dimension A[30,p]/£ for space dimension d = 3 and for representative values of p as functions 
of anisotropy parameters Qi and a2. 




FIG. 5: Behavior of the critical dimcusion A[31,p]/£ for space dimension d = 3 and for representative values of p as functions 
of anisotropy parameters Qi and a2. 



V. CONCLUSIONS 



In this paper wc have analyzed asymptotic behavior of the structure functions Sn of passively advected vector field 
with small-scale anisotropy. To this end field-theoretic renormalization group and the operator-product expansion 
have been used in a minimal-subtraction scheme of analytic renormalization. 

It is shown that the leading-order asymptotic behavior of the structure functions is determined by the isotropic 
sector of the velocity field. At the leading order in the inertial interval all the structure functions are flat, i.e., 
independent of the separation distance, with powerlike corrections (with real positive exponents) effected by the 
small-scale anisotropy. We have calculated numerically the anomalous correction exponents up to order N = 51 
to explore possible oscillatory modulation or logarithmic corrections to the leading powerlike asymptotics, but have 
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FIG. 6: Behavior of the critical dimension A[50,p]/e for space dimension d = 3 and for representative values of p as functions 
of anisotropy parameters qi and a2- 
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FIG. 7: Behavior of the critical dimension A[51,p]/e for space dimension d = 3 and for representative values of p as functions 
of anisotropy parameters qi and 02- 



found no sign of this kind of behavior: all calculated corrections have had purely powerlike behavior. Our results 
show that the exponents of the powerlike corrections tend to grow with increasing relative impact of the anisotropy. 

From the renormalization-group point of view the present model of passively advected vector field in the presence 
of strong anisotropy has turned out to be very similar to that of passively advected scalar field _13J. In particular, the 
P functions and the one-loop contributions to renormalization matrices of relevant composite operators are the same. 
Since in the published analysis of the scalar problem T^l there were some misprints, we have also presented corrected 
complete results of the calculation of the renormalization matrices. 

However, physically the two models differ significantly: instead of the anomalous powerlike growth of the structure 
functions of the scalar problem in the inertial range, in the present vector case the leading asymptotic term turned 
out to be constant with powerlike corrections due to the small-scale anisotropy. Thus, for the passively advected 
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vector field a much stronger departure from the Kolmogorov-like scaling than in the case of passively advected scalar 
is predicted. 
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